裂缝诱导双相HTI介质模型及其弹性波传播方程   总被引：1，自引：0，他引：1  

高炜  刘财  郭智奇  兰慧田  黄航 《世界地质》2014,33(4):904-933

将Biot双相介质理论与Gurevich裂缝各向异性理论相结合,建立了能够同时考虑实际裂缝性储层孔隙性和各向异性的裂缝诱导双相HTI介质模型。从本构方程、动力学方程和动力学达西定律出发,推导出了裂缝诱导双相HTI介质中弹性波传播的一阶速度-应力方程,并针对方程的刚性问题,给出了利用显式二阶时间积分法数值求解该方程时所需要满足的稳定性条件。该方程能够定量地给出双相HTI介质的波场特征与裂缝参数、背景孔隙介质参数之间的关系,描述弹性波在这种介质中的传播机理。  相似文献   

双时间步法在二维浅水方程求解中的应用          下载免费PDF全文

喻海军  黄国如  武传号 《水科学进展》2014,25(4):542-549

针对传统算法效率低的问题,将隐式双时间步法应用于求解二维浅水方程,建立了非结构网格下高效的有限体积模型。在应用双时间步法时,虚拟时间层中的定常问题采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法进行迭代求解。通过模拟计算4个典型算例以及与传统显式算法进行比较,对模型精度、效率及处理实际问题能力进行检验,分析了时间步长、内迭代次数对模型性能的影响。结果表明,双时间步法放宽了稳定性对时间步长的限制,时间步长可取到显式格式10倍以上,计算耗时减少了50%以上,模型具有良好精度与适应性,具有较好的推广应用价值。  相似文献   

关于二阶伴随模型的理论研究   总被引：9，自引：2，他引：7  

韩桂军  何柏荣  马继瑞  李冬 《海洋学报》2000,22(3):15-19

Hesse矩阵-目标函数关于控制变量的二阶偏导数形成的矩阵,在变分数据同化过程中以及敏感性分析等方面起着重要的作用;它可以通过建立数学模型的一阶和二阶伴随模型求得.以浅水方程模型为例,利用泛函的Gâteaux微分和Hilbert空间上伴随算子的概念,导出了一阶和二阶伴随模型并由此得到Hesse矩阵.改进了Zhi Wang等(1992)建立的二阶伴随模型理论.  相似文献   

利用伴随法优化非线性潮汐模型的开边界条件Ⅰ．伴随方程的建立及“孪生"数值试验   总被引：1，自引：1，他引：0  

韩桂军  马继瑞  刘克修  李冬 《海洋学报》2000,22(6)

本研究采用基于最优控制理论的伴随法,把观测资料同化到陆架海域潮汐数值模型中去,优化开边界条件,提高数值预报的精度．潮汐模型的控制方程为考虑平流项、非线性底摩擦和侧向涡动粘性项的非线性浅水方程组;采用Lagrange乘子法建立了伴随模型．研究分两部分：第一部分即本文,建立非线性浅水方程模型的伴随方程、给出目标函数的梯度,并实现“孪生”数值试验;第二部分另文给出．  相似文献   

Second kind Chebyshev operational matrix algorithm for solving differential equations of Lane–Emden type     

《New Astronomy》2013

In this paper, we present a new second kind Chebyshev (S2KC) operational matrix of derivatives. With the aid of S2KC, an algorithm is described to obtain numerical solutions of a class of linear and nonlinear Lane–Emden type singular initial value problems (IVPs). The idea of obtaining such solutions is essentially based on reducing the differential equation with its initial conditions to a system of algebraic equations. Two illustrative examples concern relevant physical problems (the Lane–Emden equations of the first and second kind) are discussed to demonstrate the validity and applicability of the suggested algorithm. Numerical results obtained are comparing favorably with the analytical known solutions.  相似文献   

A two-dimensional finite volume hydrodynamic model for coastal areas: Model development and validation     

A.H. Daoud  K.A. Rakha  A.G. Abul-azm   《Ocean Engineering》2008,35(1):150-164

A two dimensional implicit finite volume scheme for solving the shallow-water equations is developed. The effects of the Coriolis force, surface wind stress, and waves are included. A non-uniform rectilinear forward staggered grid is used with Cartesian coordinates. The time integration is performed using the Euler implicit technique. The convective flux is treated using the deferred correction method. The viscous terms are discretized using a second order central difference approximation. The SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm is used for coupling the velocity components and the water elevation gradient for the water level correction. The system of equations is solved sequentially using the Strongly Implicit Procedure (SIP). To simulate wave driven current, a phase averaged wave model is used first to simulate wave transformation and calculate radiation stresses. The performance of the developed model is validated for different sources of external forces and different combinations of boundary conditions. The validation cases include tidal circulation in a harbor and wave induced currents behind a breakwater parallel to the coastline. The model is finally applied to simulate the flow pattern in a closed artificial lagoon and along the coastline near Damietta Port located along the Northern coast of Egypt. Results of the developed model agree well with the published results for the considered cases.  相似文献   

附有限制条件的间接平差秩亏时解法初探     

秦永宽  黄声享  张书毕 《海洋测绘》2009,29(1):7-9

在附有限制条件的间接平差的应用中，若所列观测方程中实际用到的未知参数的个数少于所选参数的总数u时，就会出现法方程系数阵秩亏的问题。为了解决此问题，在理论公式推导的基础上，给出了一种新的解决方法，并结合算例进行了验证。  相似文献   

Parameter estimation for basin-scale ecosystem-linked population models of large pelagic predators: Application to skipjack tuna   总被引：2，自引：1，他引：1  

Inna Senina  John Sibert  Patrick Lehodey 《Progress in Oceanography》2008,78(4):319-335

A Spatial Ecosystem and Population Dynamic Model (SEAPODYM) is used in a data assimilation study aiming to estimate model parameters that describe dynamics of Pacific skipjack tuna population on ocean-based scale. The model based on advection–diffusion–reaction equations explicitly predicts spatial dynamics of large pelagic predators, while taking into account data on several mid-trophic level components, oceanic primary productivity and physical environment. In order to improve its quantitative ability, the model was parameterized through assimilation with commercial fisheries data, and optimization was carried out using maximum likelihood estimation approach. To address the optimization task we implemented an adjoint technique to obtain an exact, analytical evaluation of the likelihood gradient. We conducted a series of computer experiments in order to (i) determine model sensitivity with respect to variable parameters and, hence, investigate their observability; (ii) estimate observable parameters and their errors; and (iii) justify the reliability of the computed solution. Parameters describing recruitment, movement, habitat preferences, natural and fishing mortality of skipjack population were analysed and estimated. Results of the study suggest that SEAPODYM with achieved parameterization scheme can help to investigate the impact of fishing under various management scenarios, and also conduct forecasts of a given species stock and spatial dynamics in a context of environmental and climate changes.  相似文献   

Behaviour of the SMF method for the numerical integration of satellite orbits     

Pablo Martín  José M. Ferrándiz 《Celestial Mechanics and Dynamical Astronomy》1995,63(1):29-40

The SMF algorithms were recently developed by the authors as a multistep generalization of the ScheifeleG-functions one-step method. Like the last, the proposed codes integrate harmonic oscillations without truncation error and the perturbing parameter appears as a factor of that error when integrating perturbed oscillations. Therefore they seemed to be convenient for the accurate integration of orbital problems after the application of linearizing transformations, such as KS or BF. In this paper we present several numerical experiments concerning the propagation of Earth satellite orbits, that illustrate the performance of the the SMF method. In general, it provides greater accuracy than the usual standard algorithms for similar computational cost.  相似文献   

A Time-Transformed Leapfrog Scheme     

Seppo Mikkola  Sverre Aarseth 《Celestial Mechanics and Dynamical Astronomy》2002,84(4):343-354

We present a time-transformed leapfrog scheme combined with the extrapolation method to construct an integrator for orbits in N-body systems with large mass ratios. The basic idea can be used to transform any second-order differential equation into a form which may allow more efficient numerical integration. When applied to gravitating few-body systems this formulation permits extremely close two-body encounters to be considered without significant loss of accuracy. The new scheme has been implemented in a direct N-body code for simulations of super-massive binaries in galactic nuclei. In this context relativistic effects may also be included.  相似文献